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ABS  TRACT 

It  is   shown  that  the  Dahlquist  result  limiting  the  order  of 
A-stable  multistep  methods    also   applies   to  multivalue  methods. 

INTRODUCTION 

In   1956,   Dahlquist    [l]   showed  that  the  maximum  order  of  a 
strongly  stable  k-step  linear  multistep  method  is  k+1.      In   1963,  he 
showed   [2]   that  the   maximum  order  of  an  A-stable  multistep  method  is   two. 
In   1966,   the    author    [3]   showed  that   an  extension  of  multistep  methods, 
called  multivalue  methods  ,    led  to  a  larger  class   of  methods  , 
called  modified  k-step  methods  ,  such  that  the  order  of  a  stable  method 
could  be   increased  to  2k.      Since   then,   the  question  of  whether   a  multi- 
value  method  of  order  greater  than   two   could  be  A-stable  has    remained. 
The   answer  is  no  and  is   shown  by   reducing  the  problem  to  the  problem 
studied  by  Dahlquist. 

PROOF  OF  RESULT 

Consider  the  test  equation  y'    -   Ay  =  0.      Let 

a  =    [a    ,a    , .  .  .  a  ]     =    [y  ,hy '  , . .  h  y       /k!  ]     be  the   data  saved  in  the  normal 

T  T 

form  of  a  multivalue  method,   and  let  F(a.)   =   5.-,    a,  -   H  ^   a  =    a     -   ya 

where   y  =  hA.      Then,    from  equation    (11.21)    in  reference    [U],   the  multivalue 

method  is 
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a     =  Aa     _    -  I  [—  £  ]        F(Aa     .  ) 
-n  -n-1       —3a—  -n-1 


=  a^.!  -  Llii  i  -  y  «2  -]      (^-i  "  y  ^S)A^h-i 
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Aa     _ 
-n-1 


=  Sa     „ 
-n-1 

where  l_  =    [&,£,.  ..  ]T. 

If  the  method  is  A-s  table,   then   a     =  Sa__-*Oasn->-0°   for  fixed  h  whenever 

-n  —0 

Re(y)    <  0.      Hence,  A-s t ability  implies    that   all  eigenvalues   of  S   are   inside 

the  unit   circle  when  Re(y)    <  0. 

The  eigenvalues   of  S   are   the  same   as   the  eigenvalues   of 


I  -  I 


Let  A  -    EI  =    [c^,c    ,...  ,c,  ]  where    c.    is   the   i-th    column   of  A  -    £1 ,    and  let 
A£    =  v.      Then  we  want  to  study  the  solutions   of 


0  =   det 


=   det[c_  +  v 


=  dettcQ.^.Cg,...  ,ck] 


+   det[^r"^7   '  %'C2"-"%] 


1        M    0 


"  ^Uo'^-1^    •  % S*} 
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1        ^0 


v 
^{detU-Sl)   -   det[c0,  —  ,   Cg,...,^]} 
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*1  "  ^£o 


-  vydet(A-£l)    -   det[j^   ,   c.^ ,. . .  .cj} 


[p(5)   -   ya(C)]  (1) 


where  p(£)    and  a(0    are   the  weighted  sums   of  some   determinants  which 
are  polynomials   in    £   and  independent  of  y.      Now  this   is   the  same   as   the 
polynomial  that   determines   the  roots   of  the   constant   coefficient 
difference  equation   arising  when   a  linear  multistep  method  is   applied 
to  y'    =   Ay.      If  the  method  has   order  p,  then   one   eigenvalue  of  S,  namely 
one    zero  of  (1)   must    approximate   e      to   order  p,    and  this   is   the  necessary 
and  sufficient   condition   for  p    and  a  to  determine   a  p-th  order  multistep 
method.      If  the  eigenvalues   of  S   are   less   than   one   for  all  Re(u)    <  0, 
then  the   corresponding  multistep  method  is   A-stable  ,  but  this   is  not 
possible   if  p   >  2. 
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